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Abstract 

Recent data show a large difference of the genuine three-body Bose-Einstein 
correlations in S-Pb collisions and in Pb-Pb central collisions being close to 
zero in the first case and to one in the second one. These results, unexpected 
from conventional approaches, are naturally explained by the percolation of 
colour strings produced in the collisions and subsequent incoherent fragmen- 
tation of the formed clusters. 
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The strength of the genuine three-particle Bose-Einstein correlations can be measured 
by the weight factor uj introduced in 

[cm - 1} - {C2{qi2) - 1} - (C^fe) - 1} - [C2{q,i) - 1} 
UJ = , (1) 

2^{C2(gi2) - l}{C2(g23) - l}{C2(g3i) - 1} 
Here C2 and C3 are the two- and three-body correlation functions and 

2 2,2,2 /o\ 
qij=qi~qj, q = ^12 + ^23 + %i (2) 

The weight-factor uj has been experimentally studied in e"'"e~ collisions by L3 collaboration, 
with the result consistent with u; = 1 p[. It has also been studied in heavy- ion collisions. 
NA44 collaboration ||,|] has obtained uj = 0.20 ±0.02 ±0.19 for SPb collisions, i.e. compat- 
ible with no genuine three-body correlations. On the other hand, the same experiment with 
the same methodology has found 00 = 0.85 ± 0.02 ± 0.21 for central Pb-Pb collisions. This 
value is compatible with u = 0.606 ± 0.005 ± 0.178 earlier reported by WA98 collaboration 
i- 

It has been recently suggested that an explanation of these data seems to be in 
line with the behaviour of the chaoticity parameter A which measures the strength of the 
two-body Bose-Einstein correlations 

A = C2(g = 0)-l (3) 

The intercept A is the most poorly determined parameter due to Coulomb correction, de- 
pendence on the shape of correlation functions fits and other uncertainties but experimental 
data in heavy-ion collisions show that for moderate atomic numbers of colliding nuclei A 
decreases with atomic number, as expected from the corresponding increase of the number 
of independent incoherent sources Indeed going from 0-C to 0-Cu, 0-Ag and 0-Au, 

A falls from 0.79 to 0.32 [Q. However for heavier nuclei A no longer decreases and eventually 
starts to increase. For S-Pb and Pb-Pb collisions NA44 obtains A = 0.56 and 0.59 respec- 
tively |T^. Similar values have been found at RHIC for Au-Au collisions at y/s = 130 GeV 



This behaviour can be understood assuming that in a coUision colour strings are formed 
stretched between the projectile and target, which then break due to formation of quark- 
antiquark pairs. Each colour string is assumed to have a finite transverse dimension of area 
5*1 = vrrQ (ro ~ 0.2 fm). As the energy and/or atomic number of the projectile and target 
increase the number and density of strings grows, so that they start to overlap, forming 
clusters, which act as new effective sources of particle production. Both the strings and 



their clusters can be assumed to be totally chaotic sources with A = 1 |12[. Assuming that 
for particles coming from different strings there are no Bose-Einstein correlations one 
then obtains ||14| 

A = ^ (4) 

where ns and ut are the average numbers of particle pairs produced in a given rapidity and 
transverse momentum range from the same cluster and from all the clusters respectively. 

It is clear from (^) that A decreases with the number of incoherent sources (clusters). 
For very large energies and/or atomic numbers the clusterization process will diminish the 
number of independent sources (asymptotically to unity). As a consequence the chaoticity 
parameter A will grow. 

This approach may be realized in the framework of different scenarios depending on the 
assumed form of the interaction between strings at close distances [l^JT^. If the area of 
a cluster is formed by the geometrical sum of overlapping strings then a phase transition 



is observed when the string density rj reaches a critical value rjc |jT^. This percolation 
phase transition corresponds to the appearance of at least one cluster which spans the whole 
interaction transverse area. The value of rjc lies in the interval 1.17-1.5 depending on the form 
of the profile functions of the colliding nuclei. Staying within this percolation scenario, the 
dynamics of the string interaction still admits different possibilities. The observed behaviour 
of A favours considering each cluster as a single string with a higher colour given by the 
vectorial sum of colours of the overlappig strings times a factor which takes into account the 
degree of overlapping. As a result, the number of particles /i„ produced by a cluster of area 



Sn formed by n strings is given by 



nS 

where fii is the number of particles produced by a simple string. In the case of total 
overlapping Sn = Si and /in = y/nfii. In the opposite case when strings just touch each 
other Sn = nSi and one gets = n/ii as expected. 
To calculate uj we have to know 

A3 = C3(g = 0)-l. (6) 

Under the above mentioned assumptions, similarly to (|^), we have 

A3 = 5^ (7) 

where now n'g and n'rp are the average numbers of particle triplets produced in a given 
rapidity and transverse momentum range from the same cluster and from all the clusters 
respectively. 

A completely chaotic cluster of n strings will produce (l/2)/i^ pairs of particles and 
(l/6)/i^ triplets of particles with /i„ given by (||). Summing over all formed clusters i = 
1,2, ...M one obtains ns and n'^. The total number of particles produced from all clusters 
is obviously 

M 

/^ = II/^n, (8) 

1=1 

The total numbers of pairs and triplets are (l/2)/x^ and (l/6)/i^ respectively. Thus we find 

<EZ,n^SnJS,>^ ^^_^ <j:rU{n.SnJS,f'> 
< {Efii ^Jn.SnJS,y> ' < {Efii ^JnAJS,y> 

To calculate (^) a Monte-Carlo simulation was performed. We generated discs of 
radius tq inside a circle of radius R corresponding to the interaction area. For sufficiently 
large number of discs (strings) and large R both A and A3 result depending only on the 
string density rj determined by 



Identifying all the formed clusters and determining their areas we obtained A and uj as 
functions of rj. Our results are presented in Figs. and ^. The two experimental points 
of Fig. ^ correspond to central S-Pb and Pb-Pb collisions at SPS energies. The number of 
strings for these collisions is obtained from a Monte-Carlo code ]19|] based on the Quark- 



Gluon string model with tq ~ 0.2 fm and R corresponding to the required centrality. 

As to the chaoticity parameter A, the obtained behaviour for it is similar to our previous 
calculations in |]18|. Note that there a more elaborate approach was chosen in which we 
took into account energy-momentum conservation. The energy-momentum of each string 
was determined from that of the partons at its ends, which in its turn was given by the 
corresponding structure functions. Energy-momentum conservation limits the number of 
formed strings with energy sufficient to produce particles. As a result, it shifts the minimum 
of the curve of Fig. |l| to the right and makes its rise somewhat slowlier. 

The dependence oi uj on r] is found to be stronger than for A, which was to be expected 
due to stronger dependence on the number of sources. This explains why the values of uj 
are measured to be so different for S-Pb and Pb-Pb collisions at the same energy, and this 
is the reason of the good agreement obtained. 

We expect that inclusion of energy-momentum conservation will only slightly modify the 
shape of the dependence of uj on r/, as is the case of A. Experimental information on uj at 
lower values of rj (e.g. from light nuclei or peripheral heavy ion collisions) would be most 
welcome to verify the predicted change of sign of uj. 
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FIG. 1. The chaoticity parameter A as a function of rj from Eq. (|9|). The experimental points 
are for semi-central S-Pb collisions [3] (filled triangle), 18% central Pb-Pb collisions [4] (nonfilled 
box) and 10% central Pb-Pb colhsions [20] (filled box) at SPS. 
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FIG. 2. The weight factor w as a function of rj from Eqs. (|l|) an (^). The experimental points 
are for semi-central S-Pb collisions [3] (filled triangle) and 9% central Pb-Pb collisions [4] (nonfilled 
circle) at SPS. 
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